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We consider the role of the third dimension in the conductivity of a quasi 2D 
electron gas. If the transverse correlation radius of the scattering potential is smaller 
than the width of the channel, i.e. the width of the transverse electron density 
distribution, then scattering to higher levels of the confinement potential becomes 
important, which causes a broadening of the current flow profile. The resulting 
conductivity is larger than that obtained from a 2D Boltzmann equation. A magnetic 
field, parallel to the driving electric field, effectively competes with the confining 
potential and, in the limit of a strong magnetic field, it is the field, which largely 
shapes the electron and current profile, rather than the potential. As a result the 
current flow profile increases and a negative longitudinal magnetoresistivity of the 
quasi 2D electron gas may be observed. 

I. INTRODUCTION 

Impurity-limited conductivity of a quasi-two-dimensional electron gas (quasi-2DEG) is 
usually calculated by means of a quasi-classical 2D Boltzmann equation (see, e.g. Jl|, 0, [3|. [|. 
H HI) and references therein). Quantum corrections to the 2D conductivity are assumed to 
be only due to the weak localization or interaction mechanisms. The starting point of the 
quasi-classical approach is the Hamiltonian 

+ E ^,/ 3 (q 2 )e 4q2r2 a^ k2+q2/2 a /3 , k2 _ q2/2 , (1) 

a,/3,k 2 ,q2 

where 

M Qi/3 (q 2 ) = J dz%(z)U(q 2 ,z)y p (z). (2) 
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Here U(q2,z) = J d 2 r2e~ tCi2r2 U(r2, z) where U(r2,z) is a random scattering potential; the 
subscript 2 denotes here and below the 2D in-plane vectors. E a is the a-th eigenenergy of 
the quantum well, V(z), which confines the motion of the electrons in the z direction. ty a (z) 
is the corresponding eigenf unction. 

In order to derive a quasi-classical Boltzmann equation one assumes that only the highest 
occupied level, ao, of the potential V(z) is of importance and all transitions, both real and 
virtual, to lower or higher levels can be discarded. If X z characterizes the width of the 
wave function ty ao (z), then the characteristic energy interval between the state E ao and the 

h 2 

neighboring states can be estimated as AE a « — . Then the above assumption for 

2mXj 

the real transitions is justified if the channel is narrow enough and the temperature is low 
enough, AE ao > k B T '. 

Virtual transitions to other states may be caused by the random scattering potential 
U(q2,z). Its off-diagonal matrix elements (fj) are negligible if the scattering potential is 
smooth, i.e. it is characterized by a large correlation radius, r c ^> X z . Then the scattering is 
responsible only for an in-plane relaxation of the electron momentum of the non-equilibrium 
2DEG. One may then discard any possible renormalization of the current due to an ad- 
mixture of states with a ^ ao and retain only the diagonal matrix elements M Q0)Q0 (q) in 
Eq. (p. Then a quasi-classical 2D Boltzmann equation follows straightforwardly. If one is 
interested in the profile of the current flow density along the z coordinate, it coincides in 
this case with the electron density profile, determined by the wave function of the «o level, 
i.e., j(z) ~ Pao(z) = \ip ao (z)\ 2 . The width of this profile may be called width of the channel. 

A completely different situation takes place in the opposite limit, r c <C X z , to be called 
below the quantum limit. The scattering potential [/(q2,z) induces strong transitions to 
other levels a ^ a of the quantum well. The matrix elements M Q , ji a(q 2 ) with a,P ^ a 
cannot be discarded, and the conventional approach based on a quasi-classical 2D Boltzmann 
equation is not applicable. An admixture of other states with the wave functions, localized 
in much wider regions than that of ip ao (z), may lead to a much broader z profile of the 
current density. Its decay along the z axis is characterized by a length b, which we call 
effective width of the channel. It may essentially exceed the width X z . 

The conductivity of a quasi-2DEG may be now sensitive to an external in-plane magnetic 
field, which may influence the X z value and lead to a corresponding increase of the effective 
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width of the channel, which carries the current flow. As a result, a negative magnetoresistiv- 
ity in an in-plane magnetic field is expected. It is worthwhile to distinguish this mechanism 
of the longitudinal magnetoresistivity from the other mechanism recently proposed in Ref. 
0, which considers a 2D system without any account of the third dimension. The magnetic 
field polarizes electron spins and causes a change of the Fermi energy and, hence, of the 
scattering time. It is emphasized that the mechanism, we propose here, does not consider 
electron spins or their polarization at all. It is of a crucial importance below that the system 
is quasi two-dimensional rather than really two-dimensional. 



II. 2D - CONDUCTIVITY 



A calculation of the conductivity of a quasi-2DEG in the quantum limit A 2 3> r c cannot be 
carried out within the framework of a conventional quasiclassical 2D Boltzmann equation. A 
quantum approach is necessary, which takes into account the renormalization of all relevant 
quantities due to an admixture of various levels, a ^ cxq, of the quantum well. We shall 
see below that it is equivalent to explicitly accounting for the z dependence of the current 
density flow. This analysis can be best carried out using a quantum kinetic equation for the 
Wigner function. A detailed discussion of a gauge invariant derivation and analysis of such an 
equation is presented in our papers ||, EL HH • A discussion of quantum kinetic equations can 



be found in books |TT|, 0, which present also introductions to the diagrammatic technique, 



proposed originally by Keldysh Jl3] 



We consider here a model, in which calculations can be carried out analytically. The 

2 2 
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confinement potential is harmonic, V(z) = — - — with the frequency u; m is the electron 
mass. We simplify the problem by assuming that only the ground state of the harmonic 

oscillator is occupied, «o — 0, then \ z = u = \ ; u is the amplitude of the zero- 

V muj 

point oscillations. We assume also that electrons can move freely in the xy-plane and are 
scattered by a potential, which is on the average homogeneous and isotropic in the whole 
three dimensional space. Its fluctuations are characterized by a certain correlation radius 
r c . This radius can be neglected as compared to any scales important for the motion in the 
xy plane. However, following the above discussion we shall keep a finite value of r c when 
considering the z motion of the electrons. Then 

(U(r 2 , z)U{r' 2 , z')) = W6(r 2 - r' 2 )l(z - z'). (3) 
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Here U 2 is the mean square fluctuations of the scattering potential. If the scattering po- 
tential is created by short-range defects with a concentration c, then U 2 = cU% where Ud is 
the scattering potential of a defect. The correlation function for the fluctuations in the z 
direction is chosen in the Gaussian form 

l(z - z>) = exp[-(z - zT/2r 2 c }. (4) 

Calculating various space conditional moments of the quantum kinetic equation one ob- 
tains an infinite set of equations, providing the so called hydrodynamic formulation of the 
problem (for details see fTOfl). We restrict ourselves by the s- scattering in the xy-plane, 
which implies that the equations for the kinetic moment are decoupled from the equations 
for higher powers of the electron momentum. In order to receive these three equations we 
first take Eq. (31) of our paper ||10|| , multiply it by the 3D electron momentum p, and 
integrate its left and right hand parts over the variables e and p. The resulting equation 
for the z component of the nonequilibrium part of the conditional moment is trivial with 
zeros in both left and right hand sides. It corresponds to the absence of a current in the z 
direction. The equations for the x and y components read 



eE 2 n 2 p(Z) = ^j^J T^Pa exp {teftEa [P^ - ^P] } I, (5) 

where 

n 2P (Z) = -2z J -0^G<(p,Z) (6) 

is the electron density, whose Z distribution is determined by a function p{Z) normalized 
to unity, meaning that n 2 is the in-plane electron density; Z and p are Wigner variables for 
the gauge invariant Green functions 0, [TCfl ; 

dp z 



exp|^^^-^^| 



2nh 

Tr (a z V(Z,p z )G(Z,p z ) - G(Z,p z )X(Z,p z )a z ) . (7) 

Here the Groenewold [14] notations are used, according which left and right arrows over the 

d 

differential operators = — denote the operators, which act either on the left or the right 
functions in the product, respectively. G(Z,p z ) and H(Z,p z ) are matrices of the Keldysh 
Green functions and mass operators. To shorten the notations their dependencies on the 
variables p 2 and e were suppressed. 
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The mass operator accounting for the s (in the xy plane) scattering by the potential @ 
reads 

U 2 



U f 

>( e;Z,p z ) = — — ^ J d 2 p 2 dp' z G <> (p 2 ,e;Z,p' z )l(Pz-p , z )- (8) 



where l(p z ) is the Fourier transform of the correlation function (|4]). Eq.(|5]) is solved in the 
linear, with respect to the driving electric field E2, approximation. Then the s scattering 
means that the mass operators (f|) are at equilibrium. 

It is important now to distinguish the Z and p z dependencies of the quantities in Eq.(§). 
We consider the case when the Fermi level lies in the lowest subband of the confinement 
potential V(Z), then at equilibrium 

p( z ) = Po( z ) = —7= ex P (--2) • 



G^o^e-Z^) =Gf <> (p 2 ,e)p (Z,p z \u), (9) 
where G^ <:> (p 2 , £) is the equilibrium two dimensional Green functions, and 

{u 2 p 2 I 
^2 £ f^o(^) ( 10 ) 

(see discussion of the shape of the Wigner quasi-distribution function of a harmonic oscillator 
in 0). 

The time r 2 of the electron scattering in the xy plane at equilibrium does not depend on 
the variables Z or p 2 and is defined as 

^ / ^Po(Z,p>Mp z -p' z ) =i(W<(Z, Pz ,e)-rt°>(Z,p z ,e)). (11) 
t 2 J litn 

with 1(0) = l/(r c V27r). 

From the principle of the detailed balance for the quantum collision integral |J, we 
conclude that the nonequilibrium parts of the two above Green functions coincide, i.e. 
5G > (Z, p z ; p 2 , e) = 5G < (Z,p z ;p 2 ,e) = 5G(Z,p z ;p 2 ,e). However, their Z and p z depen- 
dencies do not necessarily coincide with those of the equilibrium Green functions ([]). These 
dependencies can be found from the integral equation, which follows from Eq. (||). A Gaus- 
sian dependence of the nonequlibrium part of the Green function solves this equation. We 
look for a solution in the form 



5G(p 2 ,e; Z,p z ) = 5G 2 (p 2 ,e)p (Z,p z \b) 



(12) 
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where p (Z,p z \b) is described by Eq. (|10|), in which the length u is substituted for an 
unknown length b. Here 5G2{p2,z) is the (yet unknown) nonequilibrium correction to the 
2D Green functions, Eq. (|9|). 

We now introduce all the above assumptions into the collision integral (0), and after some 
tedious but straightforward calculations arrive at the expression 



r r „„ , , f „o / 1 r 2 



M-^zZ^^-^-wU (13) 



where 



and 



1 1 1 

u 2 u 2 b 2 



1111 

f? = rf + tf + S?- (M) 



As a result the kinetic equation (El) acquires the form 



eE 2n2?0 (Z)=2/A/||_ p2/ ^ ) . (15) 



Now we compare the Z dependence of the right hand side of the kinetic equation (fH) ) . 
determined by Eq. fll"3|), and that of the left hand side, determined by the distribution 
~p (Z). The requirement that these dependencies coincide, leads to the value 



1 1/11 



which determines the scale of the Z dependence of the nonequilibrium Green function (|12|) . 

b is now the effective width of the channel which does not necessarily coincide with the width 

u of the electron density profile. One can readily see that b w u only in the limit of the large 

u 2 

correlation radius r c ^> u. However, it may become very large, b ~ — ^> w, in the quantum 

r c 

limit, r c u. 

Eq. ([15|) allows one to find directly the first conditional moment of the nonequilibrium 
correction to the 2D-Green function 

72 



de f d p 2 N ,en 2 r 2 b 



JG ' s iwiw p * !fe,£H '^ E! (17) 

and, hence, the density of the electric current in the quasi-2DEG becomes (see Eq. (45) in 

H) 

3{z) = - 2ie Stnl4^ SG ^ e) ^ m = ff2 ^ exp {-!?)*>■ (18) 
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where a 2 = is the 2D-Drude conductivity. 

m 

Integrating Eq . fll8|) over Z one gets the conductivity of the quasi-2DEG, a = a 2 n „ , 

V2r c 

which strongly depends on the ratio of the width of the channel u and the correlation 

length r c of the scattering potential. The Drude formula is applicable, if only long range 

fluctuations are characteristic of the scattering potential, i.e., r c ^> u, then — > 1, and 

V2r c 

a = a 2 . However, in the quantum limit, r c <C u, when the scattering potential fluctuates 

in the range smaller than the width of the channel, a strong deviation from the Drude 

formula is expected. The effective width of the channel, where the current actually flows, 
u 2 

— 3> u, may become much larger than the width u of the electron density profile. As 



u 2 



a result, the conductivity may also become much larger a = cr 2 — 3> 02 



III. APPLICABILITY OF THE LINEAR RESPONSE APPROXIMATION 

The above results are obtained within the linear response approximation. In order to 
probe the applicability limits of this approximation, one considers the Z profile of the drift 
velocity. It can be found using the conditional moment ( |17D (see Eq. (43) in [|i~0|j ) 



1 r de f d P2 p 



Vd "^ = MZ)J^hJ (2#i 5G2(P2 ' £) - (19) 
This definition corresponds also to the conventional equation J(Z) = en (2')vrf T .jj t (Z), in 
which n (Z) = n 2 p (Z) is the electron density at equilibrium. 

One can readily see that differing Z dependencies of the current and the electron densities 
lead to a growth of the drift velocity with increasing Z. This fact indicates that at however 
small electric field E 2 , there are large enough distances Z, at which conditions can be 
achieved, when the linear approximation in E 2 is violated. The detailed analysis requires 
extremely cumbersome and lengthy calculations, which cannot be presented here. However, 
the most important conclusions can be made in a simpler way. We obtain a rough but 
sufficient (upper limit) criterion of the applicability of the linear approximation by requiring 
that the additional energy acquired by an electron due to the current flow does not exceed 
the temperature, i.e. VdriftPF ^ ksT where pp is the Fermi momentum. Hence, the linear 
approximation is violated if Z > Z* where 

er 2 p F u E 2 



Z = U2 2 ln 

er — u 



m r c ksT 



s 



It is important to emphasize that although this condition is obtained for the harmonic 
potential well, considered here, its meaning is more general. For any potential well exists a 
value Z* so that at larger distances the linear approximation does not hold. At Z > Z* one 
has to consider a nonlinear problem involving contributions of various inelastic scattering 
mechanism (say, electron-phonon interaction). However, without making detailed calcula- 
tions one can understand that in this region the scattering becomes much stronger and the 
current flow profile is actually much lower, than that obtained in the linear approximation, 
hence, its contribution to the total current can be neglected. When integrating over Z 
one may just cut the integration at \Z\ = Z* which provides reasonable estimates for the 
nonlinear corrections to the 2D conductivity, 




Z*/b 

Z* increases with the decreasing electric field, and at Z* ^> b these corrections are hardly 
observable. 



IV. INFLUENCE OF AN IN-PLANE MAGNETIC FIELD 



It is worthwhile to discuss here the influence of an in-plane magnetic field B on the 
conductivity. According to the intuition based on the classical ideas about the electron 
motion, we do not expect any influence of an in-plane magnetic field, especially if it is directed 
parallel to the electric field. Nevertheless, as we demonstrate below, such an influence exists 
and a positive longitudinal magnetoresistivity is expected. We have to start again from Eq. 
(31) of our paper |IJ and carry out the same procedure as described above. First, one can 
readily check that Eq. (|5|) retains its form in the case of B||E 2 and the magnetic field B does 
not appear explicitly. However, it influences the shape of the equilibrium Wigner function 
TTJ , which becomes now 



Po(Z,p z \u, v) = Airhu' 2 exp < ^ } exp 



where 



v= 1 + 



2 r wvp x 

n v J u v 

4\ 1/2 



B 



(20) 



/ 1 
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he 

and l 2 B = — is the magnetic length. We use this function in order to repeat the calculations 

en 

carried above for the case of B = 0. The solution 5G(p2,s; Z,p z ) of the quantum kinetic 
equation will be again looked in the form (p~|) where the function (|20|) with b substituted 
for u is used instead of po(Z,p z \b). 

Then we get an equation determining the parameter b as a function of the magnetic field 



u 2 z/ 3 z/ 4 — 1 u 2 



1 u 2 \ 2z/ 3 



b 2 {v) 2z/ 3 - 1 I 2z/ 3 ~ r 2 \ \ 2u 3 

As a result the density of the magnetic field dependent electric current fll8|) in the quasi- 
2DEG becomes 

3(Z, v) = a 2 * exp f-T^rl E 2 . (21) 

where 

u 2 u 2 1 . / u 2 \ 

-7+o(^ + ^ 3 ) 1 + TT • 



f2(z/) r2 2 V V &2 , 

Integrating Eq.(pT|) over Z the magnetic field dependent conductivity of the quasi-2DEG is, 

( ) = b M 




FIG. 1: Relative 2D resistivity as a function of the longitudinal magnetic field in units u 2 /lg for 
different ratios of the parameter u/r c from to 2. The curves nearly coincide. 
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Figure 1 shows the relative variation of the resistivity p[y) / p(0) = er(0) / a(u) as a function 
of the magnetic field in units u 2 jl\ for four values of the parameter u/r c : 0, 0.5, 1, 2. The 
curves for u/r c = and 0.5 are indistinguishable in this plot. Two other curves for the 
higher values of the parameter u/r c also nearly collapse. In the high magnetic field limit 
the resistivity rapidly decreases. It happens when the magnetic length Ib becomes much 
smaller than the width of the channel. Then the dynamics of electrons is controlled by the 
magnetic field and the confining potential well becomes actually irrelevant. 

Experimentally this limit can hardly be achieved. If we take a typical width of the channel 
to be 100A then the high magnetic field limit is reached only in the field measured in hundreds 
or even thousands of Tesla. It is much more realistic to expect some measurements carried 
out in the low field limit (Ib > u) where a negative magnetoresistivity can be expected. We 
can write the low field expansion of the relative magnetoresisitvity as 



p(0) Ij 



where 

2 + 7v 2 + 5v A + Vl + v 4 - v A (-A + 5y/l + v 2 ) 



k 



4V1 + v A (l + VI + v A )(-v 2 + VI + v A ) 
The dependence of this coefficient on the parameter v 2 = u 2 /r 2 is nearly linear as shown 
in Figure 2 and starts from k(0) = 3/8 at v — 0. All the above data indicate that the 
resistivity may increase by a factor of 2 if the magnetic length Ib becomes comparable with 
the width of the channel u. 

A positive longitudinal magnetoresistivity was observed in SiGe layers [|T^, The 
authors attributed the effect to the electron spin polarization (a theory is presented in 
0). This explanation seems to be quite reasonable since the electron concentrations in 
these experiments are so low that a complete spin polarization is possible. The mechanism 
proposed in this paper may become more important at higher electron concentration. We 
do not know currently about any relevant experimental data and the question remains open. 

The results of this paper emphasize that the difference between the 2DEG and quasi- 
2DEG, however subtle it is, may be of an utter importance. If the scattering potential 
comprises fluctuations with a small enough correlation radius, the profile of the electron 
current may become broader than the profile of the electron density. A longitudinal negative 
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2 4 6 8 ~ 10 



FIG. 2: The dependence of the coefficient k in the low field expansion of the magnetoresistivity 
on the parameter v 2 . 

magnetoresistivity is proposed as an experimentally observable consequence. 
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